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COHPARIbON  OF  RESULTS  OF  ATMOSPHERIC  POLLUTION 
MEASUREMENTS  WITH  STANDARDS  OF  AIR  QUALITY 

Andrzej  Kasprzyckl 

In  the  article  there  is  suggested  a method  of  quantitative 
comparison  of  the  results  of  measurements  of  the  concentration 
of  atmospheric  pollutants  with  air  quality,  for  an  arbitrary  Interval 
of  time.  The  quantitative  criteria  were  determined  on  the  basis 
of  mathematical  statistics  by  the  method  of  confidence  Intervals. 

There  were  considered  two  basic  variants;  when  the  obtained 
results  produced  a simple  test  and  when  they  produce  a stationary 
time  series.  The  problem  of  comparing  emission  measurement 
results  with  air  quality  standards  was  solved  on  the  basis  of  an 
empirical  model  worked  out  in  the  USA  by  Larsen.  The  essence 
of  the  operation  is  the  establishing  of  a guaranteed  decision 
function  for  the  two  basic  cases. 

In  the  checking  of  the  sanitary  state  of  the  atmosphere 
especially  in  urban  agglomerations  and  industrial  centers  a highly 
significant  role  is  played  by  the  referring  of  the  measurements 
obtained  of  pollutant  emissions  to  legally  establish  standards 
of  air  quality. 

In  the  case  when  the  exceeding  of  the  compulsory,  for  the 
given  area,  maximum  permissable  concentrations  of  pollutants 
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has  a consequence  a certain  action,  for  instance  in  the  form  of  ' 

a limitation  of  the  emission  of  pollutants  to  the  atmosphere,  ; 

it  is  indispensible  to  establish  quantitative  criteria  on  the  ; 

basis  of  which  there  will  be  untertaken  the  decision  about  such 
limitation,  with  a desired  level  of  guarantee. 


i In  the  continuation  there  will  be  discussed  problems  connected 

I'  with  the  establishment  of  a guaranteed  (at  a set  level  of  confl- 

p dence)  decision  function  in  the  case  where  the  standards  of  air 

I quality  are  overrun  and  where  they  are  not  overrun  at  the  given 

[ ■ place  and  time.  Similarly  the  consideration  of  two  variants 

is  indispensible  namely  where  the  obtained  measurement  results 
produce  a simple  test  where  they  produce  a time  series  of  definite 
I radius  of  correlation. 

Processing  of  Emission  Measurements  that 
^ Produce  a Simple  Test 

1 I 

In  the  case  of  making  checking  measurements  of  atmospheric 
j pollution  in  a given  area  in  sufficiently  large  time  steps 

(Intervals),  greater  than  the  radius  of  correlation  of  the  time 
series  — among  the  successive  measurements,  or  in  the  case  of 
taking  mobile  measurements  with  a continuous  change  of  the  place 
of  the  Instrument,  the  collection  of  data  obtained  can  be  handled 
as  a simple  test,  l.e.  accepting  that  there  is  no  autocorrelation 
in  the  so-called  time  series. 

' Simple  Test  of  Overrunning  of  the  Established 

' Air  Quality  Standards 

It  is  established  that  emission  measurements  are  done  with 
a time  of  gatheringC?)  T^  and  it  is  desired  to  refer  them  to  air 
quality  standards  for  the  same  Interval  of  time.  Let  further 
the  number  (size)  of  the  simple  test  (non-correlated  time  series) 
amount  to  N (and)  according  to  assumption  n (n<N)  measurements 


FTL-ID(RS)I-1708-76  2 


exceed  the  established  air  quality  standard  of  the  given  chemical 
substance . 

In  the  example  considered,  a decision  function  guaranteed  at 
an  assumed  level  of  confidence,  is  equivalent  to  the  lower  estimate 
from  the  random  test  of  the  probability  of  exceeding  the  air 
quality  standard.  Since  the  inequality  j/ > Y„  occurs  with  an 
empirical  probability  p = -,  if  then  there  is  satisfied  the 
relation  (See[l : §4 . 4]  ) . 


N 


N + t{b) 
Hh) 


(p  + 


P) 


N 


+ 


tHb) 

2N 

4N*/ 


>K 


(1) 


where 

t(b)  — is  the  abscissa  of  Student’s  distribution  for  N-1  degrees 
of  freedom  and  confidence  level  1—b , 

K — is  the  established  in  the  given  problem,  critical  value 
of  the  probability, 

y — random  variable,  in  the  problem  being  considered,  the 
arbitrary  concentration  of  the  atmospheric  pollutant. 

Then  on  the  given  protected  area  there  is  undertaken  an  action 
condltlonned  on  the  overrunning  of  the  air  quality  standard  for 
a time  of  gathering  confirmed  at  the  confidence  level  1-b . 

At  present  in  Poland  there  is  in  effect  a manner  of  interpre- 
tation of  the  legally  established  maximum  permlssable  concentrations 
of  atmospheric  pollutants,  by  some  chemical  substances,  which 
permits  the  overrunning  for  Instance  of  medium  24  hour(?)  standards 
in  the  course  of  5%  of  a day  per  year  (K=0.5),  of  a 20  minute 
standard  in  the  course  of  0.5%  of  a day  per  year  (K=0.05). 

Obviously  in  this  case  the  test  (sample)  must  be  uniform  in 
relation  to  the  year's  period,  which  in  general  is  not  satisfied. 
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In  the  established  measurement  place,  for  concentrations 
for  Instance  of  20  minutes.  It  Is  best  to  use  a so-called  random 
test  proportionately  stratified  [2].  In  such  a test  the  beginning 
element  Is  selected  randomly  from  the  first  2^  hours  for  measure- 
ment, the  next  from  the  2nd  24  hours  but  now  not  randomly,  it  is 
shifted  byaa  unit  (that  is  by  the  time  of  gathering  of  a unit 
observation  for  Instance  Tj^=20  minutes,  in  relation  to  the  time 
of  the  element  already  selected  etc.  The  test  can  also  be 
concentrated  by  taking  for  instance  k aments  at  uniform  distances 
from  each  other  from  the  first  24  hour  period  and  in  the  next 
24  hour  period  all  elements  are  shifted  by  a unit  etc. 

Such  a manner  of  interpretation  of  the  standards  in  effect 
for  air  quality,  besides  the  incomprehensibility  for  the  average 
user  of  them  (l.e.)  for  one  not  making  use  of  mathematical 
statistics)  in  the  opinion  of  the  author  does  not  secure  to  a 
sufficient  degree  the  hygienic  state  of  the  air  in  urban  agglomera- 
tions since  under  critical  meteorological  conditions,  so-called 
stagnations  (lasting  generally  for  several  days),  the  appearance 
of  high  concentrations  of  atmospheric  pollutants  is  not  sporadic 
(random)  but  constitutes  the  rule  generally.  On  the  other  hand 
it  can  also  happen  that  the  sum  of  the  periods  of  stagnation 
in  a particular  place  will  be  for  Instance  less  than  5%  of  a 
day  per  year,  and  then  in  the  theory  of  the  regulations  in  force, 
there  are  no  cases  of  overrunning  of  the  standards  in  force  of 
atmospheric  pollution.  These  regulations  moreover  Ignore  levels 
of  atmospheric  hazard  in  stagnation  periods  (each  value  of 
emission  greater  than  the  normative  quantile  of  order  for  instance 
5%  is  treated  the  same  way).  For  information  it  is  suggested 
that  for  instance  in  the  USA  there  are  standardized  only  maximum 
values  (i.  e.  extreme  values  of  the  random  test)  for  different 
times  of  gathering  of  a unit  measurement,  which  (max  vallles) 
can  be  overrun  at  the  most  once  in  the  course  of  a year  [3]* 


tiimple  Test  of  not  Overrunning  of  the  Established 
Air  Quality  Standards 

A situation  of  this  type  appears  when  on  the  basis  of  emission 
measurements  made  It  is  desired  to  have  a guarantee  (at  the  established 
confidence  level  1-b ) that  at  the  given  place  and  time  the  air  quality 
standards  are  not  overrun. 

In  the  situation  analyzed  the  conclusion  has  to  do  with  the 
quantile  which  lies  above  the  extreme  element  of  the  test;  it  is 
indispensable  then  to  make  an  assumption  which  determines  the  type  of 
statistical  distribution  of  the  population.  For  measurements  of 
atmospheric  pollution  concentration  it  is  possible  to  accept  that 
they  produce  a log-normal  distribution,  l.e.,  it  is  postulated  that 
the  transformed  random  variable: 

z = log  (y  + C) 

has  a normal  distribution.  In  the  first  approximation  it  is  possible 
to  assume  that  C = 0 and  on  the  basis  of  a random  test  (of  size  N) 
to  determine  the  parameters  of  the  distribution  from  the  equations; 


m = ^ log  <y,y2 . - . yN) 

N 

1 

log^yj  — W-m* 


(2) 

(3) 


In  the  continuation  there  is  indlspenslble  the  establishment  the  order 
K of  the  quantile  y(K),  which  in  the  analyzed  statistical  distribution 
can  not  be  greater  than  the  established  air  quality  standard  of  the 
same  gathering  time  as  the  measurements  made,  l.e.,  T^. 

In  the  case  of  the  assumption  that  the  air  quality  standards  in 
the  contemplated  measurement  period  cannot  be  in  general  overrun,  we 
can  accept  (that): 


where : 


(^) 


T — is  the  period  of  time  in  the  course  of  which  the  random  test  \ 

i 

(sample)  is  obtained.  1 


I 

When  moveover  there  is  permitted  a certain  established  percent  \ 

i 

of  time  of  overrunning  the  standard  maximum  concentrations  of  j 

atmospheric  pollutants,  as  already  was  explained  above:  K = 0.05,  ■ 

K = 0 . 5 etc . I 

i 

1 

On  the  basis  of  the  assumptions  made  it  is  obvious  (that):  i 

i 

z(K)  = log  t/(K)  i Z„  = log  Y„  (5) 

•j 

] 

where : 

— is  the  standard  value  of  atmospheric  pollution  (for  a time  of 
gathering  of  a unit  measurement  T^).  ] 

i 

1 

Further  in  agreement  with  [1:  §^.7]  with  sufficient  accuracy  for  < 

the  problem  under  consideration,  it  is  possible  to  estimate  the  asympto- 
tically normal  dispersion  of  the  distribution  of  the  quantile  z(k)>  ; 

of  order  K from  the  equation: 

i 

i 

S{K)  = /r+  0,5  t*(K)  (6) 

yiv 

where : 

N — is  the  size  of  the  sample  (test) 
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s — Is  the  magnitude  determined  by  equation  (3), 
t(h)  - Is  the  abscissa  of  the  standard  normal  distribution  for  a 
value  of  Its  ordinate  K, 

as  well  as  In  accordance  with  relatively  simple  considerations 
(specially  discussed  In  [2])  If  there  occurs: 

Z„>  s t(K)  + m + s(K)  t(a)  ( 7 ) 


where : 

m and  s — are  the  magnitudes  determined  by  the  equations  (2)  and  (3) 
t(a)  — Is  the  abscissa  of  Student's  distribution  for  N-1  degrees 
of  freedom  and  a condldence  level  1-a. 

then  at  the  confidence  level: 

l-b  = (l  — a)  (1— K)  (8) 

one  can  accept  (undertake  the  decision)  that  the  air  quality  standards 
are  not  overrun  (In  the  sense  of  the  determination  (4)  of  another  for 
Instance  K = 0.05;  K = 0.5  etc.). 

It  Is  necessary  to  note  that  equation  (8)  Is  a consequence  of  the 
assumption  of  Independence  of  the  determination  (4)  (or  for  instance 
K = 0.05  etc.)  and  of  the  estimate  of  the  quantile  of  the  distribution 
(5)  by  the  method  of  confidence  Intervals  (with  the  dispersion  (6)  and 
the  confidence  level  1-a).  Such  an  assumption  appears  natural  since 
for  Instance  in  the  case  of  knowledge  of  the  precise  values  of  the 
parameters  m and  s.  In  place  of  the  Inequality  (7)  It  Is  possible  to 
use  the  relationship:  Z„  > s t(K)  + m,  reaching  Z with  a probability  less 
than  1-K.  Since  the  confidence  Interval  1-a  Is  independent  of  the 
assumed  course  of  the  argumentation  and  It  can  be  established  arbi- 
trarily, then  if  it  Is  accepted  (that): 
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1 — b 
1 — K 


(9) 


then  it  will  be  possible  to  eliminate  the  auxiliary  magnitude  a and 
assure  the  interpretation  of  the  results  with  the  guarantee  desired 
in  the  given  problem,  i.e.,  at  the  confidence  level  1-b . 

Processing  the  Results  of  Emissions  which 
Produce  a Correlated  Time  Series 

In  many  cases  the  acceptance  of  the  assumption  about  the  mutual 
independence  of  the  results  of  the  measurements  of  pollutants  emissions 
leads  to  errors  which  are  to  large.  In  a case  of  this  kind  one  has 
to  deal  with  the  interpretation  of  mean  2H  hour  concentrations  of  pollu- 
tants obtained  in  a continuous  manner  in  the  course  of  a year  as  also 
with  the  processing  of  the  records  of  automatic  instruments. 

For  stationary  time  series  (i.e.,  which  produce  a dependent 

test  (sample)  it  is  best  to  replace  such  a series  with  an  equivalent 

to  it,  in  the  sense  of  the  amount  of  information,  simple  test  (sample) 

with  a hypothetical  size  N . For  this  purpose  there  is  indispensable 

z 

an  estimate  of  the  autocorrelation  matrix  of  the  series  determined  by 
the  relationship: 


1 r,  rj 

. . Tfc 

ri  1 n 

. . Tk-l 

»■*  ’•fc-i 

. . 1 

(10) 


In  the  first  approximation  it  is  possible  to  accept  that  the  ef- 
fective level  (order)  of  the  matrix  is  significantly  less  than  the 
number  (size)  of  the  series  (i.e.,  k<^N)  as  well  as  to  estimate  the 

elements  of  the  matrix  with  the  help  of  point  estimators  from  the 
general  relationship 
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(11) 


N — VO  N 

2 2 z,Zj 

1 = 1 J = I + to 


where : 

w = 1,  2,  . . . k 

m and  s are  determined  by  the  equations  (2)  and  (2) 

As  was  shown  In  the  works  [1]  and  [2]  the  substitute  number  (size) 
for  a simple  test  equivalent  to  the  given  series  Is  obtained  from  the 
equation : 


N 


9 fc-l 

1 + ^ 2 (fc  — D’-i 
k 1-1 


(12) 


The  whole  continuation  of  the  argumentation  Is  the  same  as  for 
simple  tests,  whereby  equations  (1)  and  (7)  remain  proper  In  the 
continuation  In  the  case  of  (by)  replacing  the  number  (size)  of  the 
actual  series  N by  the  substitute  number  (size  determined  by 
equation  (12). 

Comparison  of  Pollution  Emission  Measurement 
Results  with  Air  Quality  Standards  of  Another 
Gathering  Time 

The  concept  of  an  air  quality  standard  In  protected  areas  refers 
generally  to  atmospheric  air  which  Is  used  by  human  groups  (with  the 
exception  of  factories)  for  breathing.  Here  the  standard  maximum 
possible  concentrations  of  harmful  substances  of  relatively  short 
time  of  gathering  of  a single  measurement  (20  min.  to  2^  hours)  must 
above  all  protect  against  the  appearance  of  marked  (sudden)  Instances 
of  falling  111,  moreover  In  order  to  Insure  against  Instances  of 
becoming  chronically  111,  that  Is  due  to  atmospheric  pollution,  the 
establishment  of  standards,  monthly,  seasonal,  and  yearly,  - Is 
Indlspenslble . 
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Tlie  results  of  the  view  presented  Is  that  there  exists  a need 
for  comparison  of  atmospheric  pollutant  emission  measurement  results, 
carried  out  In  the  region  of  their  Influence  on  people  and  plants, 
with  standards  of  air  quality  with  varied  tlme-of-gatherlng  of  a 
single  measurement. 

In  the  case  of  a fixed  atmospheric  pollutant  measuring  network, 
supplied  with  continuous  analyzers,  for  example  of  the  coulometrlc 
type,  the  most  rational  manner  of  processing  the  measurement  data  for 
the  above  depicted  goal.  Is  the  using  of  a yearly  random  sample, 
proportionately  stratified.  If  such  a sample  establishes  (constitutes) 
a population  compounded  of  elements  of  gathering  time  for  Instance 

equal  20  min.  It  Is  Indlspenslble  to  establish  the  relationships  that 
make  possible  the  determination  of  the  population  parameters  (time 
series)  of  this  same  physical  phenomenon . (1 . e . , atmospheric  pollutions 
by  a given  substance)  but  of  other  times  of  gathering  T2  of  a single 
measurement,  for  Instance  one  hour,  2^  hours  etc.  For  this  purpose 
It  Is  possible  to  use  an  empirical  model  of  the  distribution  of  the 
time  Intervals  of  the  emission  measurements  suggested  by  Larsen  [3 
and  4]  and  based  from  the  theoretical  side  on  the  treatment  [2]. 

As  earlier  explained  In  the  first  approximation  Is  Is  possible 
to  assume  that  the  measurements  of  atmospheric  pollutants  on  protected 
areas,  of  arbitrary  time  of  gathering  of  a( single)  measurement  — 
produce  log-normal  statistical  distributions  (this  means  that  the 
random  variable  z = log  y has  a normal  distribution).  Keeping  In 
mind  that  among  the  geometric  mean  m and  the  mean  geometric  devla- 
tlon  Sg  as  well  as  the  parameters  of  the  distribution  m and  s, 
determined  by  equ  .tlons  (2)  and  (3)  there  Is  reached  (there  occurs) 
the  relation 


on  the  basis  of  equations  specially  worked  out  in  works  [2]  and  [3] 
there  is  obtained: 


= In 


Sj  = Si 


In  T/Tj 

w = 

In  T/r, 


where : 


and  m2  — are  average  values  of  the  logarithmic  distributions  of 
the  concentrations  of  atmospheric  pollutants  determined  by  equation 
(2)  with  time  of  gathering  of  a single  measurement  and  T2 , 

and  S2  — are  dispersions  (stan  devs?)  of  the  logarithmic  distri- 
butions of  the  atmospheric  pollutant  concentrations  determined  by 
equation  (3)  of  time  of  gathering  of  a single  measurement  and 
m — is  the  average  value  of  the  distribution  of  atmospheric  pol- 
lutants (of  an  arbitrary  time  of  gathering  of  a single  measurement 
and  not  logarithmic), 

T — base  reference  period  of  Larsen's  time  model,  regularly  one  year 
(8760  hours ) . 

In  order  to  determine  the  decision  function  guaranteed  at  the 
confidence  level  1-6,  it  is  established  that  the  hypothesis  about  the 
log-normal  distribution  of  the  pollutant  emissions,  as  also  the 
functions/relations  (15),  (l6)and  (17)  are  satisfied  precisely.  In 
such  a case  using  the  method  of  confidence  Intervals  as  well  as  the 
relationship  (6)  for  a standard  of  air  quality  of  time  of  gathering 
(accumulation)  T2  there  is  then  obtained: 

— in  the  case  when  at  tiie  level  of  confidence  1-b  the  standards  are 
overrun 


V jv ; r * 


log  Y„(T2)  < Sjf(K)  +7712  — 

,/l  + 0,5  W ■ (18) 


I 

I 

— in  the  case  when  at  the  level  of  confidence  1-b  the  standards  are 
• not  overrun: 

log  Y„{Ti)  > S2t{K)  +7712+- 

+ J/1  + (19) 

where : 

N — is  the  number  (size)  of  the  original  time  series  (of 
time  of  gathering  of  a single  observation  T^), 

and  s^  — are  the  magnitudes  determined  by  the  relations  (15) 
and  (16)  for  the  point  estimators  (2)  and  (3)  of  a time  of  gathering 
of  a single  observation  1 = 1,  2 . . . N. 

" ^n^'^2^  ~ standard's  maximum  concentrations  of  atmo- 

spheric pollutants,  average  for  the  time  interval  T^, 
t(a)  and  t(K)  — are  quantities  of  the  same  significance  as  in  (7) 
and  ( 6 ) . 

Conclusions 

On  the  basis  of  the  presented  considerations  It  is  possible  to 
formulate  the  following  more  Important  affirmations,  which  can  have 
practical  applications  in  the  control  of  the  sanitary  state  of  the 
atmosphere  in  large  urban  agglomerations. 

1.  before  approaching  the  given  protected  area  for  active  influ- 
encing of  the  magnitude  of  pollutant  emissions  to  the  atmosphere,  as 
also  during  the  process  itself  of  establishing  permissable  values  of 
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emissions,  it  is  indispensible  to  introduce  quantitative  criteria, 
on  the  basis  of  which  tactical  decisions  will  be  undertaken  (this 
means  short  term  limitations  of  the  emissions  by  some  industrial 
plants,  implementation  of  the  process  of  desulphurization  of  chimney 
gases  in  electric  generating  plants,  etc.).  Such  criteria  are  also 
indispensible  for  determining  of  the  joint  maximum  amount  of  days  in 
the  year  in  the  course  of  which  the  reduction  of  pollutants  emitted 
to  the  atmosphere  will  be  necessary.  It  seems  that  in  the  case  of 
accumulation  of  a suitably  oroad  and  representative  measurement 
material,  of  the  concentrations  appearing  on  the  given  area  of  atmo- 
spheric pollutants  (emissions),  the  most  rational  way  of  comparing 
the  obtained  results  with  air  quality  standards  is  the  method  of 
confidence  intervals,  which  introduces  as  a consequence  the  establish- 
ment of  quantitative  criteria,  expressed  by  the  equations  (1),  (7), 
(18),  and  (19). 

2.  In  the  case  of  ordering  (giving  instructions  for)  a rather 
long  time  series  (this  means  a chronologically  obtained  random 
sample)  it  is  best  to  use  a proportionally  stratified  random  sample 

as  well  as  to  use  quantitative  criteria  (1)  in  the  case  of  overrunning 
of  the  established  air  quality  standards  and  (7)  in  the  case  where 
the  standard  has  not  been  overrun. 

3.  In  case  of  necessity  of  maximum  exploitation  of  the  informa- 
tion contained  in  the  time  series,  it  is  indispensible  to  estimate 
the  order  (number,  size)  of  the  substituted  series  (determined  by 
equation  (12)),  and  right  after  to  use  the  quantitative  criteria  (1) 
and  (7).  Such  calculations  can  get  effectively  done  only  with  the 
support  of  electronic  calculating  technology  which  is  today  in 
universal  use. 

In  the  control  (checking)  of  the  hygienic  state  of  the  atmo- 
sphere an  Important  matter  is  the  comparison  of  the  obtained  measure- 
ment results  (emission  measurements)  with  air  quality  tandards  of 
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various  gathering  times.  In  applications  there  is  recommended  the 


experimentally  tested  model  of  Larsen  as  well  as  the  quantitative 
criteria  (obtained  by  the  method  of  confidence  intervals)  expressed 
by  the  relationship  (l8)  in  case  where  the  air  quality  standards 
(of  a different  unit  measurement  gathering  time  than  the  original 
time  series)  are  overrun,  as  well  as  by  (19)  in  the  case  where  the 
standards  are  not  overrun. 
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